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ABSTRACT 

By means of fully kinetic simulations, we investigate electron acceleration 
during magnetic reconnection in a nonrelativistic proton-electron plasma with 
conditions similar to solar corona and flares. We demonstrate that reconnection 
leads to a nonthermally dominated electron acceleration with a power-law energy 
distribution in the nonrelativistic low-/? regime but not in the high-/? regime, 
where /? is the ratio of the plasma thermal pressure and the magnetic pressure. 

The accelerated electrons contain most of the dissipated magnetic energy in the 
low-/? regime. A guiding-center current description is used to reveal the role of 
electron drift motions during the bulk nonthermal energization. We hnd that the 
main acceleration mechanism is a Fermi-type acceleration accomplished by the 
particle curvature drift motion along the electric held induced by the reconnection 
outhows. Although the acceleration mechanism is similar for different plasma /?, 
low-/? reconnection drives fast acceleration on Alfvenic timescales and develops 
power laws out of thermal distribution. The nonthermally dominated acceler¬ 
ation resulting from magnetic reconnection in low-/? plasma may have strong 
implications for the highly efficient electron acceleration in solar hares and other 
astrophysical systems. 

Subject headings: acceleration of particles — magnetic reconnection — Sun: 
hares — Sun: corona 
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Introduction 


Magnetic reconnection is a fnndamental plasma process dnring which the magnetic 
held restrnctnres itself and converts its energy into plasma kinetic energies (e.g. Priest & 


Forbes (2000)). It occnrs nbiqnitonsly in laboratory, space, and astrophysical magnetized 


plasmas. An important nnsolved problem is the acceleration of nonthermal particles in the 
reconnection region. Magnetic reconnection has been snggested as a primary mechanism 


for accelerating nonthermal particles in solar hares (Masnda et al. 1994 Krncker et ah 


2010 Lin 2011), Earth’s magnetosphere (0ieroset et al. 2002 Fn et al. 2011 Hnang et al. 


2012), the sawtooth crash of tokamaks (Savrnkhin [2001 ), and high-energy astrophysical 


systems (Colgate et al. 2001 Zhang & Yan 2011). In particniar, observations of solar hares 


have revealed an efficient particle energization with 10% —50% of magnetic energy converted 
into energetic electrons and ions (Lin & Hndson 1976). The energetic particles nsnally 


develop a power-law energy distribntion that contains energy on the same order of the 


dissipated magnetic energy (Krncker et al. 2010 Oka et al. 2015). Some observations hnd 


that the emission has no distingnishable thermal component, indicating that most of the 


electrons are accelerated to nonthermal energies (Krncker et al. 2010; Krncker & Battaglia 


2014). This efficient prodnction of energetic particles poses a challenge to cnrrent theories 


of particle acceleration. 

Particle acceleration associated with reconnection has been stndied in reconnection- 


driven tnrbnlence (Miller et al. 1996 

), at shocks in the ontflow region |Tsnneta & Naito 

1998 

Gno & Giacalone 

2012), and in the reconnection layer (|Drake et al. 2006 Fn et al. 

2006 

Oka et al. 2010 

Kowal et al. 

2012 Gno et al. 2014; Zank et al. 2014). Previons 


kinetic simnlations have examined varions acceleration mechanisms dnring reconnection, 
inclnding the Fermi-type mechanism in magnetic islands (hnx ropes in three-dimensional 
simnlations; Drake et al. (2006); Gno et al. (2014[)) and direct acceleration in the dihnsion 


region (Pritchett 2006 Hnang et al. 2010). Most simnlations focns on regimes with plasma 


/3 > 0.1 with no obvions power-law distribntions observed. It was argned that particle loss 


from the simnlation domain is important for developing a power-law distribntion (Drake 


et al. 2010). Early simnlations of relativistic reconnection showed power-law distribntions 


within the X-region (Zenitani & Hoshino 2001). Recent kinetic simnlations with a highly 
magnetized (a = ^j ^ 1) pair plasma fonnd global power-law distribntions withont 
the particle loss, althongh the loss mechanism may be important in determining the spectral 


index (Gno et al. 2014, 2015). It is nnknown whether or not this is valid for reconnection in 


a nonrelativistic proton-electron plasma, since its property is different from the relativistic 


reconnection (Lin et al. 2015). 


Motivated by the resnlts of relativistic reconnection, here, we report fnlly kinetic simn- 
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lations of magnetic reconnection in a nonrelativistic proton-electron plasma with a range of 
electron and ion /3e = A = 0.007—0.2. The low-/? regime was previously relatively unexplored 
due to various numerical challenges. We hnd that reconnection in the low-/? regime drives ef- 
hcient energy conversion and accelerates electrons into a power-law distribution f{E) ~ E~^. 
At the end of the low-/? cases, more than half of the electrons in number and 90% in en¬ 
ergy are in the nonthermal electron population. This strong energy conversion and particle 
acceleration led to a post-reconnection region with the kinetic energy of energetic particles 
comparable to magnetic energy. Since most electrons are magnetized in the low-/? plasma, 
we use a guiding-center drift description to demonstrate that the main acceleration process 
is a Fermz-type mechanism through the particle curvature drift motion along the electric 
held induced by fast plasma hows. The development of power-law distributions is consistent 
with the analytical model ( Guo et ah]|2014 ). The nonthermally dominated energization may 
help explain the efficient electron acceleration in the low-/? plasma environments, such as 
solar hares and other astrophysical reconnection sites. 


In Section we describe the numerical simulations. In Section we present simulation 
results and discuss the conditions for the development of power-law distributions. We discuss 
and conclude the results in Section IH 


2. Numerical simulations 


The kinetic simulations are carried out using the VPIC code (Bowers et ah 2008), 


which solves Maxwell’s equations and follows particles in a fully relativistic manner. The 
initial condition is a force-free current sheet with a magnetic held B = i?o tanh(z/A)a; -|- 
Bq sech.{z/X)y, where A = A is the half thickness of the layer. Here, A is the ion inertial 
length. The plasma consists of protons and electrons with a mass ratio mil me = 25. The 
initial distributions for both electrons and protons are Maxwellian with uniform density 
uq and temperature fcTj = kTe = O.OlmgA. A drift velocity for electrons Ue is added to 
represent the current density that satishes the Ampere’s law. The initial electron and ion 
/?e = A = STTriokTe/B^ are varied by changing Upe/^ce, where ujpe = \//me is the 
electron plasma frequency and Gee = ^BQjimeC) is the electron gyrofrequency. Quantities 
A = 0.007, 0.02, 0.06, and 0.2 correspond to ojpe/Vtee = 0.6, 1, a/S and \/T0, respectively. 
The domain sizes are L^x = 200A x 100A- We use N^. x = 4096 x 2048 cells with 200 
particles per species per cell. The boundary conditions are periodic along the x-direction, 
perfectly conducting boundaries for helds and rehecting boundaries for particles along the 


2 r-direction. A long wavelength perturbation is added to induce reconnection (Birn et ah 


2001 ). 
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Simulation results 


Under the influence of the initial perturbation, the current sheet quickly thins down to 
a thickness of ~ de (electron inertial length c/upe) that is unstable to the secondary tearing 
instability ( Daughton et ah|2009 ; Liu et ah|2013b ). Fig.[^ a) and (b) show the evolution of the 
out-of-plane current density. The reconnection layer breaks and generates a chain of magnetic 
islands that interact and coalesce with each other. The largest island eventually grows 
comparable to the system size and the reconnection saturates at thlci ~ 800. Fig. [^c) shows 
the time evolution of the magnetic energy in the x-direction (the reconnecting component) 
and the kinetic energy of electrons and ions Ki for the case with = 0.02, respectively. 
Throughout the simulation, 40% of the initial Sbx is converted into plasma kinetic energy. Of 
the converted energy, 38% goes into electrons and 62% goes into ions. We have carried out 
simulations with larger domains (not shown) to confirm that the energy conversion is still 
efficient and weakly depends on system size. Since the free magnetic energy overwhelms the 
initial kinetic energy, particles in the reconnection region are strongly energized. Eventually, 
Ke and Ki are 5.8 and 9.4 times their initial values, respectively. Fig. [^d) shows the ratio of 
the electron energy gain AKe to the initial electron energy Ke{0) for different cases. While 
the /3e = 0.2 case shows only mild energization, cases with lower give stronger energization 
as the free energy increases. 

The energy conversion drives strong nonthermal electron acceleration. Fig. [^a) shows 
the final electron energy spectra over the whole simulation domain for the four cases. 
More electrons are accelerated to high energies for lower-/? cases, similar to earlier simu¬ 


lations (Bessho & Bhattacharjee 2010). More interestingly, in the cases with /?e = 0.02 and 
0.007, the energy spectra develop a power-law-like tail f{E) ~ E~'p with the spectral index 
p ~ 1. This is similar to results from relativistic reconnection (Guo et al. 2014, 2015). We 


have carried out one simulation with rrii/me = 100 and /3e = 0.02 and find a similar electron 
spectrum. In contrast, the case with = 0.2 does not show any obvious power-law tail, con¬ 


sistent with earlier simulations (Drake et al. 2010). The nonthermal population dominates 


the distribution in the low-/? cases. For example, when we subtract the thermal population 
by fitting the low-energy distribution as Maxwellian, the nonthermal tail in the /?e = 0.02 
case contains 55% of the electrons and 92% of the total electron energy. The power-law tail 
breaks at energy Eb ~ lOEth for /?e = 0.02, and extends to a higher energy for /?e = 0.007. 
Fig. I^b) shows the fraction of nonthermal electrons for different cases. For /?e = 0.007, the 
nonthermal fraction goes up to 66%, but it decreases to 17% for /?e = 0.2. Fig. [^c) and (d) 
show nacc/ue at tVLci = 125 and 400 for the case with /?e = 0.02, where Uacc is the number 
density of accelerated electrons with energies larger than three times their initial thermal 
energy, and Tie is the total electron number density. The fraction of energetic electrons is 
over 40% and up to 80% inside the magnetic islands and reconnection exhausts, indicating a 
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Fig. 1.— Out-of-plane current density for the case with /3e = 0.02 at (a) tflci = 62.5 and (b) 
tflci = 400. (c) The energy evolution for the (3e = 0.02 case. Ebxit) is the magnetic energy of 
the reconnecting component. Se is the electric energy. Ki and K^. are ion and electron kinetic 
energies normalized by e:bx(0), respectively, (d) The ratio of electron energy gain AiFg to the 
initial Ke for different (3e- 
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bulk energization for most of electrons in the reconnection layer. The energetic electrons will 
eventually be trapped inside the largest magnetic island. The nonthermally dominated dis¬ 
tribution contains most of the converted magnetic energy, indicating that energy conversion 
and particle acceleration are intimately related. 

To study the energy conversion, Fig. [^a) shows the energy conversion rate dSc/dt from 
the magnetic held to electrons through directions parallel and perpendicular to the local 
magnetic held. We dehne dSc/dt = ■ EdV, where V indicates the simulation domain 

and j' is j’n or We hnd that energy conversion from the perpendicular directions gives 
~ 90% of the electron energy gain. By tracking the trajectories (not shown) of a large 
number of accelerated electrons, we hnd various acceleration processes in the dihusion region. 


magnetic pile-up region, contracting islands, and island coalescence regions (Hoshino et ah 


2001 

Hoshino 2005 

Drake et ah 

2006 

Dahlin et ah 

2014 

Guo et ah 

2014 

)• 


back and forth through a Fermi-like process accomplished by particle drift motions within 
magnetic structures (Li et ah 2015, in preparation). To reveal the role of particle drift 
motions, we use a guiding-center drift description to study the electron energization for the 
(3e = 0.02 case. The initial low (3 guarantees that this is a good approximation since the 
typical electron gyroradius pe is smaller than the spatial scale of the held variation (~ di). 

By ensemble averaging the particle gyromotion and drift motions, the perpendicular 


current density for a single species can be expressed as (Parker 1957 Blandford et ah 2014) 

Bx{B \/)B 


j± = P\\ 




P. (F) X VB 


V X 


P±B 

F2 


E X B B duE 

+ P ^ I- PmJ^ X 


F2 


dt 


using a gyrotropic pressure tensor V = P±I + (Fy — P±)bb, where P\\ = rUe J fvjdv and 
Fj_ = 0.5me / fvj_dv, p is the particle charge density, and pm is the particle mass density. 
The terms on the right are due to curvature drift, VB drift, magnetization, E x B drift, 
and polarization drift, respectively. The expression is simplihed as + jg + + 

Jexb P jpj which j exb direct contribution to the energy conversion. This gives 

an accurate description for j ^ if the pressure tensor is gyrotropic. To conhrm this, we 
calculate the electron pressure agyrotropy AQe = 2 , where Fj_ei and P±e 2 are the 

two pressure eigenvalues associated with eigenvectors perpendicular to the mean magnetic 


held direction (Scudder & Daughton 2008). AQe measures the departure of the pressure 


tensor from cylindrical symmetry about the local magnetic held. It is zero when the local 
particle distribution is gyrotropic. Fig. |^b) shows that the regions with nonzero A0e are 
localized to X-points. The small A0e indicates that the electron distributions are nearly 
gyrotropic in most regions. Therefore, the drift description is a good approximation for 
electrons in our simulations even without an external guide held, which is required for this 
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Fig. 2.— (a) Electron energy spectra f{E) at tQci = 800 for different /3e- The electron 
energy E is normalized to the initial thermal energy Et^. The black dashed line is the 
initial thermal distribution, (b) Time evolution of the fraction of nonthermal electrons for 
different initial /3e. rinth is fhe number of nonthermal electrons obtained by subtracting the 
htted thermal population from the whole particle distribution. The fraction of electrons 
with energies larger than three times the initial thermal energy at (c) tQci = 125 and (d) 
= 400. 



















Fig. 3.— (a) Energy conversion rate dSc/dt for electrons throngh the parallel and perpen- 
dicnlar directions with respect to the local magnetic field, compared with the energy change 
rate of electrons dKg/dt for the case with = 0.02. The shown valnes are integrals over the 
whole simulation domain, (b) Electron pressure agyrotropy A0g at = 400 in the same 
case. See the text for details. 


description in a high-/3 plasma (Dahlin et al. 2014). 


Fig-i a) and (b) show time-dependent deddt and Eg from different current terms, where 
Eg = j^{dEc/dt)dt. The contribution from polarization current and parallel current are small 
and not shown. The curvature drift term is a globally dominant term of j^ ■ E, the VB term 
gives a net cooling, and the magnetization term is small compared to these two. Fig. |^c) 
shows the spatial distribution of jg ■ E. When the flow velocity u is along the magnetic 
held curvature n due to tension force, jg ■ E ^ {P\\B x k/B'^) ■ {—u x B) > 0. These 
regions are a few di along the 2 ;-direction, but over SOd* along the x-direction. The overall 
effect oi jg ■ E is a strong electron energization. Fig. |^d) shows that jg ■ E is negative in 
most regions because the strong Vi? is along the direction out of the reconnection exhausts. 
Then, j ■ E r\j {B X VB) ■ {—u X B) < 0. Note at some regions, j ■ E can give strong 


acceleration. Fig. |^e) shows the cumulation of the jg ■ E and jg ■ E along the x-direction. 
In the reconnection exhaust region(x = 60 — llSdj), jg ■ E is stronger than jg ■ E, so 
the electrons can be efficiently accelerated when going through these regions. In the pile- 
up region(x = 120 — IdOdj), k, VB and u are along the same direction, so both terms 
give electron energization. In the island coalescence region(x ~ ISOdj), jg ■ E gives electron 
heating, while jg ■ E gives strong electron cooling. Although the net effect is electron cooling. 


island coalescence can be efficient in accelerating electrons to the highest energies (Oka et al. 


2010 ). 


It has been shown that the curvature drift acceleration in the reconnection region cor¬ 
responds to a Fermi-type mechanism (Dahlin et al.||2014 Guo et aT]|2014 2015). To develop 
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Fig. 4.— Analysis using a drift description for the case with /3e = 0.02. (a) The energy 
conversion rate due to different types of current terms, compared with the electron energy 
change rate dKe/dt. j^ ■ E, j^ ■ E, and ■ E represent energy conversion due to curvature 
drift, VB drift, and magnetization, respectively, (b) The converted magnetic energy due to 
various terms in (a), normalized to the initial magnetic energy of the reconnecting component 
£6x(0). (c) Color-coded contours of energy conversion rate due to curvature drift at f = 
400f2“^. K and u indicate the directions of the magnetic held curvature and the bulk how 
velocity, (d) Color-coded contours of energy conversion rate due to VB drift at t = 40014“/. 
B and VB indicate the directions of the magnetic held and the gradient of |.B|. Both ■ E 
and jg ■ E are normalized to the 0.002nQmeC^u;pe- (e) The cumulation of ■ E (blue) and 
jg ■ E (green) along the a:-direction. The black line is the sum of these two. 































10 


a power-law energy distribution for the Fermi acceleration mechanism, the characteristic 
acceleration time Tacc = 1/a needs to be smaller than the particle injection time Tinj (Guo 


et ah 2014, 2015), where a = {1/e){de/dt)^ and de/dt is the energy change rate of particles. 


To estimate the ordering of acceleration rate from the single-particle drift motion, consider 
the curvature drift velocity Vc = v‘^^B x K,/{QceB) in a curved field where Rc = 
so the time for a particle to cross this region is ~ Rc/'V\\ and the electric field is mostly 
induced by the Alfvenic plasma flow E ~ —va x B/c. The energy gain in one cycle is 
6e ~ mvAV\\. The time for a particle to cross the island is Lisiand/v\\- Then, the acceleration 
rate de/dt ~ sva/ Lisiand for a nearly isotropic distribution. The characteristic accelera¬ 
tion time Tacc ~ Lisiand/vA- Taking Lisiand ~ SOdj and va ~ 0.2c, the acceleration time 
Tacc ~ 25011^^. The actual acceleration time may be longer because the outflow speed will 
decrease from va away from the X-points, and the VB term gives a non-negligible cool¬ 
ing effect. Our analysis has also found that pre-acceleration and trapping effects at the 
X-line region can lead to more efficient electron acceleration by the Fermi mechanism and 
are worthwhile to investigate further (Hoshino 2005 Egedal et ah 2015; Huang et ah 2015). 
Taking the main energy release phase as the injection time Tinj ~ 8000“/^ the estimated 
value of Tinj/Tacc 3.2, well above the threshold. For the case with /dg = 0.2, the ratio 
Tinj/Tacc ~ 0.32 < 1, SO there is no power-law energy distribution. 


4. Discussion and conclusion 


Nonthermal power-law distributions have rarely been found in previous kinetic simula¬ 
tions of nonrelativistic magnetic reconnection ( Drake et al.||2010 ). We find that two essential 
conditions are required for producing power-law electron distribution. The hrst is that the 
domain should be large enough to sustain reconnection for a sufficient duration. A power- 
law tail develops as the acceleration accumulates long enough {Tinj/Tacc > !)• The second 
condition is that plasma (3 must be low to form a nonthermally dominated power-law dis¬ 
tribution by providing enough free energy (oc l//d) for nonthermal electrons. Assuming 10% 
of magnetic energy is converted into nonthermal electrons with spectral index p = 1, one 
can estimate that /dg is about 0.02 for half of the electrons to be accelerated into a power 
law that extends to This agrees well with our simulation. We point out that a loss 


mechanism or radiation cooling can affect the final power-law index (Fermi 1949; Guo et ah 


2014) of nonthermal electrons. Gonsequently, including loss mechanisms in a large three¬ 


dimensional open system is important, for example, to explain the observed power-law index 
in solar flares and other astrophysical processes. Another factor that may influence our re¬ 
sults is the presence of an external guide held Bg. Our preliminary analysis has shown that 
the Fermi acceleration dominates when Bg < Bq. The full discussion for the cases includ- 
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ing the guide field will be reported in another publication (Li et ah 2015, in preparation). 
A potentially important issue is the three-dimensional instability, such as kink instability 
that may strongly infiuence the results. Unfortunately, the corresponding three-dimensional 
simulation is beyond the available computing resources. We note that results from three- 
dimensional simulations with pair plasmas have shown development of strong kink instability 


but appear to have no strong infiuence on particle acceleration (Guo et ah 2014; Sironi & 


Spitkovsky 2014). The growth rate of the kink instability can be much less than the tearing 


instability for a high mass ratio (Daughton 1999), and therefore the kink instability may be 
even less important for electron acceleration in a proton-electron plasma. 


In our simulations, the low-/? condition is achieved by increasing magnetic field strength 
(or equivalently decreasing density). We have carried out low-/3 simulations with the same 
magnetic field but lower temperature and found a similar power-law distribution (Li et al. 
2015, in preparation). 


The energy partition between electrons and protons shows that more magnetic energy 
is converted into protons. For simulations with a higher mass ratio rrij/me = 100, the 
energetic electrons still develop a power-law distribution, and the fraction of electron energy 
to the total plasma energy is about 33%, indicating that the energy conversion and electron 
acceleration are still efficient for higher mass ratios. Our results show that ions also develop a 
power-law energy spectrum for low-/3 cases and the curvature drift acceleration is the leading 
mechanism. However, the ion acceleration has a strong dependence on the mass ratio mj/rue 
for our relatively small simulation domain (~ lOOdj). We therefore defer the study of ion 
acceleration to a future work (Li et al., 2015, in preparation). 


The energetic electrons can generate observable X-ray emissions. As nonthermal elec¬ 
trons are mostly concentrated inside the magnetic islands, the generated hard X-ray flux can 


be strong enough to be observed during solar flares in the above-the-loop-top region (Ma- 


suda et al. 1994; Krucker et al. 2010| and the reconnection outflow region (Liu et al. 


2013a). The nonthermal electrons may also account for the X-ray flares in the accretion 


disk corona (Galeev et al. 1979 Haardt et al. 1994 Li & Miller 1997) 


In summary, we find that in a nonrelativistic low-/3 proton-electron plasma, magnetic 
reconnection is highly efficient at converting the free energy stored in a magnetic shear into 
plasma kinetic energy and accelerating electrons into nonthermal energies. The nonthermal 
electrons contain more than half of the total electrons, and their distribution resembles power- 
law energy spectra with spectral index p ~ 1 when particle loss is absent. This is in contrast 
to the high-/? case, where no obvious power-law spectrum is observed. It is important to 
emphasize that the particle acceleration discussed here is distinct from the acceleration by 


shocks, where the nonthermal population contains only about 1% of particles (Neergaard 
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